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ABSTRACT 

This paper explores the effect of rotation on the K-mechanism instability of slowly 
pulsating B stars. A new nonadiabatic code, that adopts the so-called ‘traditional ap¬ 
proximation’ to treat the Coriolis force, is used to investigate the influence exerted 
by rotation over the stability of stellar models covering the mass range 2.5 Mq ^ 
M, ^ 13.0 Mq. The principal finding is that, for all modes considered apart from the 
prograde sectoral class, rotation shifts the K-mechanism instability toward higher lu¬ 
minosities and effective temperatures; these shifts are accompanied by broadenings in 
the extent of instability strips. Such behaviour is traced to the shortening of mode pe¬ 
riods under the action of the Coriolis force. Instability strips associated with prograde 
sectoral modes behave rather differently, being shifted to marginally lower luminosities 
and effective temperatures under the influence of rotation. 

The implications of these results are discussed in the context of the observational 
scarcity of pulsation in B-type stars having significant rotation; various scenarios are 
explored to explain the apparent dichotomy between theory and observations. Fur¬ 
thermore, the possible significance of the findings to Be stars is briefly examined. 

Key words: stars: oscillation - stars: rotation - stars: variables: other - stars: early- 
type - stars: emission-line. Be - methods: numerical 


1 INTRODUCTION 

The slowly pulsating B (SPB) stars are a homogeneous class 
of mid-B type, main-sequence objects that exhibit multiperi- 
odic light and line-profile variations over timescales on the 
order of 1-5 days. Waelkens (1991) was the first to clas¬ 
sify these stars as a distinct group of early-type nonradial 
pulsators, by combining under the same aegis the photo¬ 
metric variables discovered by Waelkens & Rufener (1985), 
and the 53 Per spectroscopic variables first observed by 
Smith (1977). The theoretical pulsation characteristics of 
these objects have since been studied extensively, resulting 
in a canonical picture (see, e.g., Pamyatnykh 1999, and ref¬ 
erences therein) of high-order g-mode pulsation, driven by 
the same iron-group k mechanism responsible for the insta¬ 
bility of the P Cep pulsators (Dziembowski & Pamyatnykh 
1993). 

One outstanding issue in the understanding of SPB 
stars is the effect rotation has on their instability. Photo¬ 
metric surveys of open clusters by Balona (1994), Balona 
& Koen (1994) and Balona & Laney (1995, 1996) failed 
to find any evidence for SPB-like variability, but the au- 
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thors did note that the stars observed were characterized 
by moderate or rapid projected equatorial rotation veloci¬ 
ties (usini > lOOkms”^). This result has led some to sug¬ 
gest that rotation suppresses the K-mechanism excitation 
of g modes (see, e.g., Balona 2001). However, an alterna¬ 
tive interpretation has been explored by Townsend (2003b), 
who examined the effect of the rotation-originated Coriolis 
force on the flux perturbations produced by low-frequency 
g modes. He demonstrated that the confinement of these 
modes within an equatorial waveguide (e.g., Lee & Saio 
1990; Townsend 2003a) reduces the amplitude of the result¬ 
ing photometric variations, possibly to below the sensitivity 
of the observations. 

This paper addresses the SPB-rotation question from 
a more-theoretical perspective, by examining the effect the 
Coriolis force has on the K-mechanism instability of SPB 
stars. The high-order g modes excited in these objects are 
characterized by low frequencies, and their dynamics can 
thus be expected to be influenced appreciably by the Cori¬ 
olis force (see Townsend 2005). For this reason, canonical 
approaches that treat the rotation as a small perturbation 
to the pulsation (e.g., Carroll & Hansen 1982; Unno et al. 
1989, their §19 and references therein) are not well suited 
to SPB stars. Building on a methodology originated by Lee 
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& Saio (1987a,b), alternative treatments have emerged over 
the past decade (see, e.g., Lee 1998, 2001); however, these 
are characterized by such a high degree of mathematical and 
numerical complexity, that their large-scale application is - 
with present-day computational resources - impractical. 

The approach adopted herein aims at a compromise be¬ 
tween sophistication and practicality, by employing an ap¬ 
proximate method to treat the pulsation-rotation interac¬ 
tion (Sec. 2) that is well-suited to the g modes found in SPB 
stars. This approach is applied to a range of B-type stellar 
models (Sec. 3), with the results of these stability calcula¬ 
tions presented in Section 4. The findings of the analysis are 
then discussed and summarized in Section 5. 


2 METHOD 

2.1 Theoretical Treatment 

From a theoretical standpoint, the introduction of rotation 
significantly complicates the analysis of stellar pulsation. 
The centrifugal force tends to distort the equilibrium star 
away from a spherical configuration, and the Coriolis force 
leads to mixing between the radial and angular components 
of the fluid momentum. Both of these processes mean that 
the linearized equations describing nonradial, nonadiabatic 
pulsation are no longer separable in all coordinates, leading 
to a significant increase in the computational effort neces¬ 
sary for their solution (see, e.g., Lee 2001). 

However, for the case of the low-frequency g modes 
found in SPB stars, an approximate treatment of the Cori¬ 
olis force can be used to restore the separability of the pul¬ 
sation equations, and thereby render the problem tractable 
even with modest computing resources. At the heart of this 
approach is the so-called ‘traditional approximation’, intro¬ 
duced in a geophysical context by Eckart (1960), and first 
applied to stellar pulsation by Lee & Saio (1987a). The tra¬ 
ditional approximation involves neglecting the polar com¬ 
ponent of the rotation angular frequency vector ft, in the 
Coriolis terms of the linearized momentum equation. As dis¬ 
cussed by Lee & Saio (1997), this approximation is reason¬ 
able throughout regions where the pulsation frequency a 
(as measured in the corotating frame) and the rotation fre¬ 
quency are both very much less than the Brunt-Vaiasala 
frequency A/” associated with the gravitational stratification 
of the medium. In the case of g modes in B-type stars, this 
condition is fulfilled in almost all of the stably stratified, 
radiative envelope where the modes are propagative. 

In addition to the traditional approximation, the treat¬ 
ment presented herein requires a number of other simplifi¬ 
cations to permit the separation of the pulsation equations 
in all three spherical-polar coordinates These sim¬ 

plifications, which follow the analysis by Townsend (2003b), 
are: 

(i) The assumption that the rotation^ is slow enough for 
both the centrifugal force, and the departures from spheric¬ 
ity engendered by it, to be neglected. In this context, ‘slow’ 


^ Assumed throughout to be uniform. 


may be interpreted as the condition fl^, where 


flc 


8GM* 

27i?3 


( 1 ) 


defines the critical rotation frequency within the Roche 
model, with M* and i?* being the stellar mass and radius, 
and G the gravitational constant. 

(ii) The Cowling (1941) approximation, whereby pertur¬ 
bations to the gravitational potential are neglected. This ap¬ 
proximation is reasonable for all but low-order, low-degree 
modes. 

(hi) The nonadiabatic radial flux (NARF) approxima¬ 
tion, whereby the divergence of the horizontal Eulerian flux 
perturbation is neglected in the energy equation. As demon¬ 
strated by Townsend (2003b), this approximation - first in¬ 
troduced by Savonije et al. (1995) - is valid under the same 
conditions (a, O <C A/”) as the traditional approximation it¬ 
self. 

(iv) The steady-wave approximation, whereby the imag¬ 
inary part (Ji of the pulsation frequency a is neglected in 
solving the angular parts of the pulsation equations. This 
approximation, which is discussed by Townsend (2000), is 
appropriate when the growth rate of the pulsation is slow, 
so that Ui is very much smaller than the corresponding real 
part (Tr of the pulsation frequency. 


This collection of approximations and assumptions 
closely follows those employed by Ushomirsky & Bildsten 
(1998), with the exception that their use of the quasi- 
adiabatic approximation is replaced herein by the adoption 
of the NARF approximation. This replacement is particu¬ 
larly significant, in that the quasi-adiabatic approximation 
can lead to incorrect results when applied to B-type stars 
(see, e.g., Dziembowski et al. 1993); in contrast, the NARF 
approximation is able to capture the essential physics of the 
nonadiabatic processes responsible for the K-mechanism in¬ 
stability. 


2.2 Pulsation Equations 

Within the approximations outlined in the preceding sec¬ 
tion, the spatial and temporal dependence of perturbed vari¬ 
ables is expressed in the corotating frame as 

/'(r,M,<^,t) = /'(r)er(l^)e‘(’"'^+"‘> (2) 

and 

Sfir, n, cj), t) = <5/(r) 07(^) , (3) 

where Eulerian and Lagrangian perturbations are repre¬ 
sented by primes (') and 5, respectively; t is the temporal 
coordinate; and /r = cosO. Here, the Hough function 07(71) 
is a solution to Laplace’s tidal equation (see, e.g., Bildsten 
et al. 1996), that depends implicitly on the ‘spin parameter’ 
V = 2fl/(Tr characterizing the relative strength of Coriolis 
and buoyancy forces (Townsend 2005). In the limit of no ro¬ 
tation, 07 in) approaches the normalized Legendre function 
PR in) (see Lee & Saio 1990) having harmonic degree i and 
azimuthal order m. 

Inspired by the dimensionless formulations introduced 
by Saio & Cox (1980) and Lee & Saio (1997), which them¬ 
selves are built on the seminal Dziembowski (1971) treat¬ 
ment, the radial dependencies f'{r) and Sf{r) of perturbed 
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variables are expressed in terms of a set of eigenfunctions qi 
(i = 1, 2, 5, 6 )^, such that 


1 

H 

II 

p'(r) 

r 

pgr 

^^0 „ j 

5Ln{r) 

Crt 

L* 


(72 X 


= q6X 


e+i 


(4) 

(5) 


Here, x = r/R, is the dimensionless radial coordinate, 
the radial fluid displacement, p' the Eulerian pressure per¬ 
turbation, 5S the Lagrangian perturbation to the specific 
entropy, and SLb. the Lagrangian perturbation to the radial 
part of the radiative luminosity. Other symbols appearing 
in these expressions are defined by Unno et al. (1989), with 
the exception of the effective harmonic degree 


e = 


+ 4A7 - 1 


( 6 ) 


introduced by Townsend (2000) as the solution of the equa¬ 
tion i{l -b 1 ) = XT- In these relations, XT is the eigenvalue 
of Laplace’s tidal equation associated with the Hough eigen¬ 
function OTiu)- This eigenvalue depends implicitly on the 
spin parameter v, and therefore on both the pulsation fre¬ 
quency a and the rotation frequency fl. 

The radial eigenfunctions qi are found as the solutions 
of four coupled, first-order differential equations. 


, dgi 

da; 



1 5i + 

■|(l-bl) 

V 

fr. ) 


Fi 


(72 + VtX (75, 


(7) 


= (ciw^ — A*)q\ -b —U -b 3 -£)g 2 + vrx^q^, 


. d(75 

da; 


= H [Vad(t/ - ClCd") - 4(Vad - V) -b C 2 ] % -b 


V 


^^^^(Vad - V) - C2 


ClCU^ 


S1 + 

a;2 ^ 


( 8 ) 


( 9 ) 


I TV (4 - ns) - l\q5 - j — xqe, 


HV 


d(76 _ _ qi 

^ j — q 

dx X'^ 


SadCsH -b 


'-bl) dfa 
" ^ Ax 


ClCd^ 


(72 


+ 


(cses - iaacd) ^ -b 1 j ge- 


( 10 ) 


Here, the dimensionless pulsation frequency ui has the usual 
definition 


u) = a 



( 11 ) 


and all other symbols follow the nomenclature of Unno et al. 
(1989), with the exception of the introductions 



C3 = C3 Zr 




( 12 ) 

(13) 


^ The indices 3 and 4 are traditionally reserved for the perturba¬ 
tion to the gravitational potential and its radial derivative; due 
to the adoption of the Cowling (1941) approximation, they are 
not used in the present analysis. 


and 


C4 = C4 /r = 


Airr'^pTcp GM> 


m 


(14) 


In the limit —> 0, the effective harmonic degree i 
approaches the true harmonic degree i, and the governing 
equations (7-10) reduce to those describing nonradial, nona- 
diabatic pulsation of a non-rotating star, within the NARF 
and Cowling (1941) approximations. However, even allow¬ 
ing for the alternative nomenclature, the energy conserva¬ 
tion equation ( 10 ) in this limit appears rather different than 
in other nonadiabatic treatments (compare with, e.g., Unno 
et al. 1989, their eqn. 24.12®). In particular, the derivative 
of the dimensionless radiative luminosity Zr appears in the 
present treatment due to the choice of ‘frozen convection’ 
approximation; this choice centres around neglecting per¬ 
turbations to the convective source term in the linearized 
energy equation, as per eqn. 21.7 of Unno et al. (1989). Al¬ 
though Pesnell (1990) has argued that such an approach 
is unphysical, it is required to ensure the pulsation equa¬ 
tions remain self-consistent at the origin. In any case, since 
g modes do not penetrate far into convective regions, any 
errors introduced by this convection-freezing choice should 
be minimal. 

Solutions of the pulsation equations (7-10) are required 
to satisfy boundary conditions at the centre and at the sur¬ 
face. The inner conditions 


CiU!^ qi — tq 2 = 0 
3Zr gi — x^qe = 0 


0 , 


(15) 


ensure that solutions remain finite and regular at the origin. 
Likewise, the outer boundary conditions 


gi - 52 = 0 

(2 - 4VadU) qi + 4VadU 52 + 4 55 - ge = 0 

(16) 

follow from the assumptions that the surface pressure tends 
to zero, and that there is no flux incident from outside the 
star (see, e.g., Unno et al. 1989). Finally, the arbitrary over¬ 
all scaling of solutions is constrained by the normalization 
condition 



51 = 1 at a; = 1 . 


(17) 


2.3 Implementation 

To solve the pulsation equations (7-10) and accompanying 
boundary conditions (15-17), a wholly-new Fortran 95 code 
was developed. The code, named boojum, follows the root¬ 
finding approach pioneered by Castor (1971) and Osaki & 
Hansen (1973): one of the boundary conditions is set aside, 
allowing solution of the equations to be achieved at arbitrary 
dimensionless frequency u>. The suppressed boundary con¬ 
dition is then used to construct a discriminant D{u}), whose 
roots correspond to the eigenfrequencies of the full system 
of equations. 


® Note that the sign of the C3ICluj‘S term in their equation is 
incorrect. 
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In the present case, the outer mechanical boundary con¬ 
dition is used to form the discriminant 


D{lo) = ^ 

(qi + « 72 )„=o 


(18) 


The numerator is the boundary condition itself, while the 
denominator - which is guaranteed by the inner boundary 
conditions (15) never to be zero - ensures that the discrim¬ 
inant remains well behaved (see Townsend 2000). To calcu¬ 
late the solution vector qi required to evaluate this discrim¬ 
inant, BOOJUM uses a standard relaxation approach (e.g.. 
Press et al. 1992); because the pulsation equations are lin¬ 
ear in qi, only a single iteration is required at each value of ui. 
The f-dependent scalings, adopted in the definitions (4) of 
the eigenfunctions, successfully prevent the loss of precision 
near the centre discussed by Takata & Lofller (2004). For 
improved accuracy, centred finite differences are used in the 
relaxation algorithm; however, in regions where C4 > 10^, 
BOOJUM switches to one-sided differences for the thermal 
equations (9-10), in accordance with the Sugimoto (1970) 
prescription for avoiding numerical stability (and see also 
Unno et al. 1989, their §24). 

Solution of the characteristic equation 


D{io) = 0, (19) 

defining the eigenfrequencies of the pulsation equations, is 
accomplished in BOOJUM using the robust Traub (1964) im¬ 
plementation of the root-finding algorithm devised by Muller 
(1956); this algorithm is a generalization of the complex se¬ 
cant approach favoured by Castor (1971). Note that the lat¬ 
ter author used approximate roots of the characteristic equa¬ 
tion to find starting points for the solution, via iterative re¬ 
laxation, of the full pulsation equations. A somewhat simpler 
approach is adopted by BOOJUM, whereby equation (19) is 
solved to the desired fractional tolerance in uj (in the present 
work, 10~^° for both real and imaginary parts), and no ad¬ 
ditional calculations are performed. 

For calculation of the eigenvalue, required (c.f. 
eqn. 6) to evaluate the £ terms appearing in the pul¬ 
sation equations and boundary conditions, BOOJUM uses 
the matrix-mechanical approach pioneered by Lee & Saio 
(1987a). The implementation largely follows the procedure 
described by Townsend (2003b); however, the truncation di¬ 
mension N is determined dynamically, by repeatedly dou¬ 
bling N until the fractional change in X’P drops below some 
specihed threshold (taken to be 10“® throughout the present 
work). Furthermore, the algorithm due to Kahan (1966), as 
implemented by the lapack subroutine library (Anderson 
et al. 1999), is used for calculating the matrix eigenvalues; 
this algorithm performs significantly better at finding iso¬ 
lated eigenvalues than the previously-adopted QL approach. 

As input, BOOJUM is supplied with the desired mode 
parameters {£,m), the rotation angular frequency O, and an 
indication of the region over which to search for eigenfre¬ 
quencies satisfying the characteristic equation (19). The co¬ 
efficients appearing in the pulsation equations and boundary 
conditions are evaluated using a precomputed stellar model 
(see Sec. 3); the model - typically composed of ca. 1,000 ra¬ 
dial points - is interpolated onto a new grid having 10,000 
points, distributed non-uniformly via an approach similar to 
that described by Christensen-Dalsgaard & Mullan (1994, 
their Appendix A3). Cubic splines are used for interpolat¬ 



log ^eff (K) 


Figure 1. The theoretical HR diagram for the stellar models 
introduced in Section 3; these are plotted as points in the effec¬ 
tive temperature (Teff) versus stellar luminosity (L*) plane. The 
dashed lines running diagonally from top left to bottom right 
indicate the ZAMS and TAMS main sequence limits, while the 
asterisk shows the location of the 53 Per-like model analyzed in 
Section 4.1. Full evolutionary tracks for four selected models are 
also shown in the diagram, plotted as solid lines and labeled by 
their corresponding stellar masses, in solar units. The two higher- 
mass tracks do not reach the low-temperature boundary of the 
diagram, because the Warsaw-New Jersey evolutionary code is 
unable to follow these models beyond the ignition of core helium 
burning. 


ing all variables apart from the Brunt-Vaiasala frequency 
A/”; the latter is interpolated linearly, to avoid the introduc¬ 
tion of spurious oscillations in the molecular weight gradient 
zone situated outside the convective core. 


3 STELLAR MODELS 

The Warsaw-New Jersey stellar evolution code is used to 
calculate 115 tracks of stellar models, sampling the initial 
mass range M* = 2.5 Mq-5.2 Mq at a resolution 0.05 Mq, 
the range M* = 5.2Mq-9.2Mq at a resolution 0.1 Mq, 
and the range M* = 9.2-13.0 Mq at a resolution 0.2 Mq; 
each track extends from zero-age main sequence (ZAMS) 
to terminal-age main sequence (TAMS). Details of the code 
have already been given by Dziembowski & Pamyatnykh 
(1993) and Dziembowski et al. (1993); the only signihcant 
difference in the present work is the adoption of more-recent 
OPAL tabulations for opacity (Iglesias & Rogers 1996) and 
equation of state (Rogers et al. 1996). In all cases, the initial 
hydrogen and metal mass fractions are set at the canonical 
values A = 0.7 and Z — 0.02, with a heavy-element mixture 
taken from Grevesse & Noels (1993). No account is taken 
of the effect of rotation on the stellar evolution; although 
Pamyatnykh (1999) has demonstrated that rotation-induced 
mixing can shift the red edge of the SPB instability strip to 
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Table 1. Fundamental parameters of the 53 Per-like stellar model 
analyzed in Section 4.1. 

Tefj/K log3/cms-2 L*/Lq M,/Mq Rt/KQ 

15,300 4.04 600 4.80 3.48 

lower effective temperatures, the present analysis is more 
concerned with the dynamical influence of rotation on SPB 
pulsation. 

Each of the calculated evolutionary tracks is composed 
of approximately 42 stellar models, with a grand total of 
4,876 models for the entire set. Fig. 1 plots the positions 
of these models in the theoretical Hertzsprung-Russell (HR) 
diagram, along with the loci defining the ZAMS and TAMS 
boundaries. The dense model coverage of the main sequence 
band is to facilitate the accurate positioning of instability 
strips (Sec. 4.2) in the Tefi — T* plane. 

4 STABILITY CALCULATIONS 
4.1 53 Per Model 

From the models introduced in the preceding section, one is 
selected as having parameters close to the values inferred by 
De Ridder et al. (1999) for 53 Per, the archetypal SPB star; 
these parameters are documented in Table 1, and the cor¬ 
responding position in the HR diagram is shown in Fig. 1 
by an asterisk. This 53 Per-like model is employed to ex¬ 
amine the general influence the Coriolis force exerts over 
£ = 1 and 1 = 2 nonradial g modes. Eigenfrequencies of 
these modes are calculated using boojum, for all possible 
values —£ ^ m ^ £ of the azimuthal order, and over a range 
of angular frequencies ranging from the non-rotating case up 
to the intermediate rate fl/flc = 0.5. This upper limit cor¬ 
responds to a rotation period H = 1.26 days, with an equa¬ 
torial velocity Ve that is 35% of the star’s critical velocity 
Vc = 419kms“^; by way of comparison, the most rapidly- 
rotating SPB star known to date, HD 1976, rotates at ~ 32% 
of its critical velocity (Mathias et al. 2001). 

The results from these calculations are shown in Fig. 2, 
where the corotating frame period P = 2Ti/a-c of each mode 
is plotted as a function of rotation angular frequency. Points 
associated with eigenfunctions qi having the same radial or¬ 
der fS are linked together into a single, continuous curve; 
the weight of the curves, at each value of 17, is used to indi¬ 
cate whether a mode is stable (ni > 0) or unstable (ci < 0) 
against K-mechanism excitation. 

At 17 = 0, the eigenfrequencies are degenerate in az¬ 
imuthal order m, owing to the arbitrariness of the model’s 
polar axis. This degeneracy is lifted upon the introduction 
of rotation, with a splitting initially following the first-order 
relation 

A(Tr « ml7(7n,r, (20) 


^ As defined by Unno et al. (1989, their eqn. 17.5), within the 
generalization to the Cowling (1941) nomenclature introduced by 
Scuflaire (1974) and Osaki (1975). 


Table 2. Instability ranges for g modes of the 53 Per-like model; 
for each (£, m) pair considered, and at three differing rotation 
rates, the range of \n\ spanned by unstable modes is tabulated, 
along with the corresponding values of the corotating frame pe¬ 
riod P, the inertial frame period P; and the harmonic degree ratio 
£/£. The absolute value of h is used, because this radial order is 
negative for g modes. 


(£, m) 

1 

n\ 

P (days) 

Pi (days) 

£/£ 




Q,/Clc = 

0.00 




(L-1) 

12 

- 22 

1.10 

- 1.93 

1.10 

- 1.93 

1.00 - 

1.00 

(LO) 

12 

- 22 

1.10 

- 1.93 

1.10 

- 1.93 

1.00 - 

1.00 

(Ll) 

12 

- 22 

1.10 

- 1.93 

1.10 

- 1.93 

1.00 - 

1.00 

(2,-2) 

12 

- 29 

0.64 

- 1.46 

0.64 

- 1.46 

2.00 - 

2.00 

(2,-1) 

12 

- 29 

0.64 

- 1.46 

0.64 

- 1.46 

2.00 - 

2.00 

(2,0) 

12 

- 29 

0.64 

- 1.46 

0.64 

- 1.46 

2.00 - 

2.00 

(2,1) 

12 

- 29 

0.64 

- 1.46 

0.64 

- 1.46 

2.00 - 

2.00 

(2,2) 

12 

- 29 

0.64 

- 1.46 

0.64 

- 1.46 

2.00 - 

2.00 




fi/fic = 

0.25 




(1,-1) 

13 

- 19 

1.40 

- 2.07 

0.90 

- 1.14 

0.79 - 

0.75 

(1,0) 

12 

- 25 

1.03 

- 1.77 

1.03 

- 1.77 

1.09 - 

1.32 

(1,1) 

12 

- 31 

0.85 

- 1.40 

1.27 

- 3.16 

1.40 - 

2.26 

(2,-2) 

12 

- 28 

0.68 

- 1.58 

0.44 

- 0.70 

1.84 - 

1.74 

(2,-1) 

12 

- 30 

0.64 

- 1.40 

0.51 

- 0.90 

1.99 - 

2.18 

(2,0) 

12 

- 33 

0.61 

- 1.28 

0.61 

- 1.28 

2.13 - 

2.70 

(2,1) 

12 

- 35 

0.58 

- 1.18 

0.76 

- 2.21 

2.22 - 

3.18 

(2,2) 

12 

- 33 

0.58 

- 1.23 

1.07 

- 46.50 

2.25 - 

2.83 




fi/fic = 

0.50 




(1,-1) 

13 

- 18 

1.51 

- 2.10 

0.69 

- 0.79 

0.72 - 

0.69 

(1,0) 

12 

- 29 

0.89 

- 1.46 

0.89 

- 1.46 

1.32 - 

2.00 

(1,1) 

12 

- 38 

0.64 

- 1.04 

1.30 

- 5.76 

1.99 - 

4.02 

(2,-2) 

12 

- 27 

0.71 

- 1.59 

0.33 

- 0.45 

1.75 - 

1.66 

(2,-1) 

12 

- 33 

0.61 

- 1.25 

0.41 

- 0.63 

2.12 - 

2.77 

(2,0) 

12 

- 38 

0.54 

- 1.00 

0.54 

- 1.00 

2.46 - 

4.17 

(2,1) 

12 

- 42 

0.50 

- 0.85 

0.82 

- 2.59 

2.71 - 

5.58 

(2,2) 

12 

- 38 

0.51 

- 1.00 

2.73 

- -1.71 

2.60 - 

4.18 


or alternatively 

AP « -mHUft/, (21) 

derived by Ledoux (1951); here, H = 2tt/Q. is the rotation 
period, and A denotes the change in the indicated quantity. 
In the case of the high-order g modes considered herein, 
the term Cn,i approximates to !/£(£ -I- 1), explaining why 
the £ = 1 modes exhibit steeper gradients AP/Af7 around 
17 « 0 . 

Toward larger values of 17, departures from the above 
linear relations are increasingly apparent in Fig. 2. For 
all apart from the prograde sectoral (PS) modes having 
£ = — m, these departures produce an overall trend of de¬ 
creasing pulsation period with increasing rotation rate, ir¬ 
respective of whether the mode is prograde (m < 0), ret¬ 
rograde (m > 0) or zonal (m = 0). Such behaviour comes 
from the influence of the Coriolis force, and can readily be 
understood by recalling that - at the most general level - 
the frequency of a wave may be expressed as the square root 
of the ratio between a generalized stiffness and a generalized 
inertia (see, e.g., Lighthill 1978). With the introduction of 
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1= 1, m= - 1 



m-l 



n/fic 


t=l, m=0 



i-2, m—-2 



n/Oo 


Figure 2. Pulsation periods for £ = I and i = 2 g modes of the 53 Per-like model, plotted as a function of rotation frequency up to 
a maximum Q = 0.5 flc. Each curve corresponds to a particular radial order h; the line weight is used to indicate whether modes are 
stable (thin) or unstable (thick) against the K-mechanism instability. To improve the clarity of the figure, only radial orders |n| ^ 75 are 
plotted. 


rotation, the stiffness® usually due to buoyancy is augmented 
by the Coriolis force; this leads to an increase in frequency, 
and thence a corresponding decrease in mode period. 

Because the strength of the Coriolis force varies with 
spin parameter v, the effects described are differential: the 
periods of long-period, larger-j/ modes are shortened by a far 

® That is, the restoring force on displaced fluid elements. 


greater degree than those of short-period, smaller-!/ modes. 
Accordingly, the density of the pulsation spectrum, as mea¬ 
sured by the number of modes either per period interval 
or per frequency interval, increases markedly with rotation 
rate. As can be seen in Fig. 2, this spectral compression is 
especially pronounced for the m = 1 modes. 

This latter result warrants some explanation. For non- 
PS modes in the inertial regime |z/| > 1, the Coriolis force 
acts to confine the modes within an equatorial waveguide. 
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Figure 2 - continued 


whose boundaries are situated at /r = ±l/|z/| (see, e.g., Bild- 
sten et al. 1996). A requirement of the trapping is that the 
s + 1 nodes of the Hough function 07“ be fitted within these 
waveguide boundaries, where 

s = t-\m\±l (22) 


1990). In the limit \v\ 1, the fitting requirement is em¬ 

bodied in the angular characteristic equation 

\T t=^v'^{2s + lf (23) 

for the Hough-function eigenvalue (Townsend 2003a, his 
eqn. 38), corresponding to an effective harmonic degree 


is the meridional order introduced by Townsend (2003a). 
In this expression, the difference between the prograde and 
zonal modes (m ^ 0; minus sign) and the retrograde modes 
(m > 0; plus sign) comes about because an extra pair of 
nodes appears in the latter when > 1 (see Lee & Saio 




{21 — 2|m| — l)v (m ^ 0), 
{21 — 2\m\ +S) v {m > 0). 


(24) 


Regarding the ratio iji as a measure of the severity of equa¬ 
torial confinement, this latter expression indicates m = 1 
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modes - for each value of ^ - are the most affected by the 
Coriolis force. This result can be seen in the right-most col¬ 
umn of Table 2 introduced below; physically, it follows di¬ 
rectly from the fact that these modes have the largest merid¬ 
ional orders (s = 1 for f = 1, and s = 2 for I = 2,) and are 
therefore the most compressed in the polar direction when 
squeezed into the equatorial waveguide. 

The foregoing analysis does not apply to the PS 
modes, = (1,-1) and {£,m) = (2,-2), whose peri¬ 

ods lengthen somewhat toward larger O. As discussed by 
Townsend (2003a), this class of mode is transformed by the 
Coriolis force into equatorial Kelvin waves. Such waves have 
different properties than the usual gravito-inertial waves 
found in a rotating, stratified system; in particular, they 
are characterized by geostrophic balance, whereby the Cori¬ 
olis force arising from azimuthal fluid motions is countered 
by polar pressure gradients (see, e.g. Gill 1982). Due to this 
balance, the Coriolis force makes relatively little difference 
to the generalized stiffness of PS modes, and therefore does 
not produce the marked period decrease seen for the other 
modes. 

As a prelude to the analysis presented in the follow¬ 
ing section, the focus now turns briefly toward the stability 
of the modes plotted in Fig. 2. Once again, there is a di¬ 
chotomy between the PS modes and the others: the former 
are partly stabilized by rotation, and the latter partly desta¬ 
bilized. However, in both cases the rotation does not alter 
the property that a contiguous sequence of modes is un¬ 
stable toward the k mechanism. With this result in mind. 
Table 2 succinctly summarizes the data shown in the hg- 
ure, by indicating the range of the instability, at three dif¬ 
fering rotation rates, for each {l,m) pair considered. This 
range is expressed in terms of the radial orders of the 
highest- and lowest-frequency unstable modes; the corre¬ 
sponding corotating frame period P, inertial frame period® 
Pi = 2n I {a — mfl), and harmonic degree ratio ijl are also 
tabulated. From the table, it is evident that a grand total 
of 7 out of the 29 initially unstable PS modes are stabilized 
by the rotation. In the same way, the 94 non-PS modes that 
are initially unstable are augmented, upon the introduction 
of rotation, by the destabilization of a further 58; the largest 
gains, of 16 additional unstable modes, are accorded to the 
{e,m) = (1,1) group. 

An important result clearly seen in Fig. 2, is that the 
stabilization or destabilization of modes, with varying ro¬ 
tation rate, does not suffice to maintain the boundaries of 
the instability at a constant frequency or period. For in¬ 
stance, the {1,'rn) = (1,1) modes mentioned above exhibit 
instability over the period range ~ 1.05-1.95 days in the 
non-rotating case; but this range is narrowed and shifted 
to ~ 0.65-1.05 days in the rapid-rotation limit. Thus, even 
though more modes are unstable in this limit, they occupy 
a narrower region of the pulsation spectrum; such behaviour 
is a direct consequence of the spectral compression discussed 
previously. 

The processes selecting which modes are unstable are 
similar to those operative in a non-rotating star. As dis¬ 


® For retrograde modes, negative values of Pi - e.g., in the final 
row of Table 2 - indicate the mode is prograde in the inertial 
frame. 


cussed by Dziembowski et al. (1993), one requirement for 
the K mechanism to work efficiently is that the relative La- 
grangian pressure perturbation 5p/p be large and slowly 
varying with radius within the excitation zone. This con¬ 
dition establishes the short-period limit of the instability at 
|h| ~ 12 (c.f. Table 2); because the shape of radial eigen¬ 
functions is insensitive to the strength of the Coriolis force 
(see, e.g., Ushomirsky & Bildsten 1998), this limit does not 
vary appreciably with rotation rate. 

The corresponding long-period limit is determined not 
so much by the efficiency of the K-mechanism driving, but 
more by the onset of significant thermal damping. This 
damping arises due to radiative diffusion between neigh¬ 
bouring fluid elements, operating primarily in a region of the 
envelope (logT ~ 5.9) situated below the excitation zone. 
For the diffusion to be effective in stabilizing the pulsation 
requires a combination of long pulsation periods and high 
radial orders, the latter serving to steepen the temperature 
gradients driving the diffusion. The interplay between these 
two factors leads to a long-period instability limit that is 
quite sensitive to the rotation rate, as Fig. 2 attests. 


4.2 Instability Strips 


The focus is now broadened, from the 53 Per-like model 
considered in the preceding section, to the complete set of 
stellar models introduced in Section 3. At three differing ro¬ 
tation rates, fl/flc = 0.0, 0.25 and 0.5, boojum is used to 
search for unstable £ = 1 and £ = 2 g modes; the results of 
these calculations are presented in Fig. 3, where the same 
HR diagram shown in Fig. 1 is overplotted by the insta¬ 
bility strip associated with each {£, m) combination at each 
rotation rate. These instability strips enclose all models that 
are unstable (ui < 0) toward the excitation of one or more 
modes of the indicated type. 

As with the preceding section, there is a clear division 
between the PS and the non-PS modes. The effect of rota¬ 
tion on the former is to shift their instability strips along the 
main sequence, toward lower temperatures and luminosities. 
The converse behaviour is exhibited by the latter, for which 
instability strips are shifted toward higher temperatures and 
luminosities. These shifts affect the blue edges of the insta¬ 
bility to a greater degree than the corresponding red edges. 
For instance, the stellar models lying along the ZAMS are, 
in the absence of rotation, unstable toward (£,m) = (1,1) 
modes over the temperature range 4.06 logTefj 4.21. 
Upon the introduction of rotation, at the rate fl/Qc = 0.5, 
this range shifts and broadens to 4.08 Si logTefr 4.28. 

The behaviour of both classes of mode can be under¬ 
stood with reference to the prerequisites for excitation dis¬ 
cussed by Dziembowski et al. (1993). In addition to the 
aforementioned restriction on the pressure eigenfunction 
(c.f. Sec. 4.1), it is necessary that the thermal timescale in 
the excitation zone be comparable to, or longer than, the 
pulsation period. This latter condition is the key to the be¬ 
haviour seen in Fig. 3. The thermal timescale 


T-th(r) 


//*• TCpdMr 

T. 


(25) 


evaluated at the radius where the opacity derivative kt is 
maximal, is a strong function of stellar effective tempera¬ 
ture; it declines rapidly toward larger values of Tefr, because 



Influence of the Coriolis force on slowly pulsating B Stars 9 


1= 1, m= - 1 
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Figure 3. Instability strips in the HR diagram, for £ = 1 and £ = 2 g modes of the stellar models introduced in Fig. 1. In each panel, 
corresponding to a particular combination (£, m) of mode parameters, the extent of the K-mechanism instability is indicated using thick 
lines, at the three differing rotation rates Q/Qc = 0.0 (solid), 0.25 (dotted) and 0.5 (dashed) considered. The uneven boundaries of the 
instability strips, evident in some panels, are due to the discrete spacing of the stellar models in the Teff—L* plane. 


the excitation zone^ is situated ever closer to the surface. 
Consequently, for the non-PS modes, whose periods become 
shorter under the influence of the Coriolis force, the effect of 
rotation is to shift the instability toward stars having higher 
Teff, so that the Tth~P matching condition can still be met. 
Conversely, for the PS modes, the marginal lengthening in 

^ Typically occurring at a temperature logT ^ 5.3. 


periods due to rotation pushes the instability toward slightly 
lower Teff. 

Fig. 4 illustrates these processes, by plotting instabil¬ 
ity regions in the Teff-P plane for stellar models situated 
along the ZAMS and TAMS boundaries of the main se¬ 
quence. These regions, the complements of the instability 
strips shown in Fig. 3, indicate how the range of unstable¬ 
mode periods evolves with changing effective temperature; 
also shown in the figure, as separate curves for the ZAMS 
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Figure 3 - continued 


and the TAMS models, is the thermal timescale within the 
excitation zone. For the non-PS modes, the period shorten¬ 
ing caused by rotation is evidently responsible for the shift of 
the K-mechanism instability toward higher effective temper¬ 
atures, in order to maintain the loose correspondence be¬ 
tween Tth and P. Likewise, for the PS modes, the period 
lengthening causes the instability to shift to lower effective 
temperatures, once again to keep rth in step with P. 

It is instructive to relate the foregoing analysis to pre¬ 
vious studies (e.g., by Balona & Dziembowski 1999) of the 
effect that a varying harmonic degree £ has on the insta¬ 


bility of SPB stars. It can be recognized that the pulsation 
equations (7-10) are - modulo a number of approximations 
- no different than those describing modes of harmonic de¬ 
gree i in a non-rotating star. Accordingly, the influence of 
the Coriolis force on the period and stability of an individual 
mode can be followed simply by allowing I to vary in some 
specihed manner. For non-PS modes, this variation assumes 
the asymptotic form given in equation (24); clearly, any in¬ 
crease in V, corresponding to more-rapid rotation, leads to 
a proportional increase in £. But, as the instability strips 
presented by Balona & Dziembowski (1999, their Fig. 5) re- 
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Figure 4 - continued 


5 DISCUSSION & SUMMARY 

In the preceding sections, an approximate method for treat¬ 
ment of the Coriolis force (Sec. 2.1) is used to devise equa¬ 
tions governing nonadiabatic, nonradial pulsation of rotat¬ 
ing stars (Secs. 2.2). These equations are solved (Sec. 2.3) 
for a range of mid-B type stellar models (Sec. 3); the general 
Ending (Sec. 4) is that the Coriolis force shifts the instabil¬ 
ity strip associated with K-mechanism excitation of g modes, 
toward higher luminosities and effective temperatures for 
non-PS modes, and toward lower luminosities and effective 
temperatures for prograde sectoral modes. 

An immediate corollary of this result is that the Coriolis 
force is unable to stabilize all B-type stars against g mode 


pulsation; rather, it can only alter which particular stars 
- at a given rotation rate - are unstable. Similar conclu¬ 
sions were reached by Ushomirsky & Bildsten (1998) in their 
quasi-adiabatic analysis, and by Lee (2001) in his complete 
treatment of the Coriolis force. However, as discussed in Sec¬ 
tion 1, observations of SPB stars reveal an apparent paucity 
of objects rotating at signihcant rates. How might this dis¬ 
crepancy between theory and observations be resolved? 

One possibility is that the centrifugal force, neglected 
in the present treatment and by Lee (2001), can act to stabi¬ 
lize modes. Certainly, not taking this force into account will 
introduce a certain level of error in the results - especially 
at the upper limit Q/fl = 0.5 adopted herein, which comes 
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very close to violating the assumption (c.f. Sec. 2.1) that 
^ However, Lee & Baraffe (1995) and Lee (1998) 
have demonstrated that the centrifugal force has little ef¬ 
fect on the stability of the g modes characteristic to SPB 
stars. Perhaps a more plausible scenario is that resonant 
coupling, between modes of differing harmonic degree (see, 
e.g. Chandrasekhar & Lebovitz 1962), may act to suppress 
the K-mechanism instability. These couplings, which are ne¬ 
glected in the present analysis because of the adoption of the 
traditional approximation, were shown by Lee (2001) to in¬ 
hibit the excitation of selected modes that would otherwise 
be unstable in rotating SPB stars. 

Nevertheless, it is highly unlikely that resonant coupling 
could suppress the instability of all modes. With this in 
mind, it appears increasingly probable that, as the theoret¬ 
ical analysis indicates, g modes are excited in rotating SPB 
stars; but that the equatorial-confinement effects demon¬ 
strated by Townsend (2003b) render them very difficult to 
detect observationally. In support of this interpretation for 
the observational scarcity of rotating SPB stars, two of the 
open clusters originally thought to be devoid of SPB stars - 
NGC 4755 (Balona & Koen 1994) and NGC 6231 (Balona & 
Laney 1995) - have, upon closer scrutiny, been found to con¬ 
tain a number of candidate-SPB objects (see Stankov et al. 
2002; Arentoft et al. 2001). 

In addition to their direct applicability to SPB stars, 
the results presented herein are highly relevant to the un¬ 
derstanding of the Be phenomenon. The multiperiodic spec¬ 
troscopic and photometric variations exhibited by many Be 
stars (e.g., p Gen - Rivinius et al. 1998; to GMa - Stefl 
et al. 2003) are usually interpreted as arising from nonradial 
g modes (see Rivinius et al. 2003, and references therein). 
It is not unreasonable to suppose the excitation of these 
modes to be due to the same k mechanism operative in the 
SPB stars. However, a historical problem with this stance 
has been that many variable Be stars have early spectral 
types B0-B3, falling blueward of the high-temperature limit 
of the SPB instability strip for low-degree pulsation (Balona 
& Dziembowski 1999). Furthermore, periodic variations are 
not seen in Be stars with types later than B5 (e.g., Baade 
1989a,b, and references therein), even though the SPB in¬ 
stability strip extends all the way down to types B8 or B9. 

Both of these difficulties are resolved by allowing for the 
influence of the Goriolis force, which - as the present anal¬ 
ysis demonstrates - has the effect of shifting K-mechanism 
instability toward earlier spectral types. Of course, this ra¬ 
tionalization can only work if the modes seen in Be stars 
are not prograde sectoral; however, in light of the analysis 
presented by Rivinius et al. (2003), who found at least 16 
of the 27 Be stars in their sample to pulsate in retrograde 
(i, m) = (2, 2) modes, this latter restriction does not appear 
to be problematic. 
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